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1. Introduction 

Two-particle correlations in AT 7 , Ad were studied extensively at RHIC and LHC. They allow 
to disentangle different sources of correlations: jets, flow, resonance decays, quantum statistics ef¬ 
fects, conservation laws, etc. Here we report preliminary results on two-particle At] A(j) correlations 
in inelastic p+p interactions at SPS beam momenta (20, 31, 40, 80, and 158 GeV/c). 


2. ArjA0 correlations 

Correlations are calculated as a function of the difference in pseudo-rapidity ( 77 ) and azimuthal 
angle ( 0 ) between two particles in the same event. 

Afl = Iflt-ml A0 = |0!-02| 

The uncorrected (raw) correlation function is calculated as: 


C raw (Arj, A0) 


<3 S(AT7,A0) 

Ndata M (Ar? , A0 ) 


(2.1) 


where 


A(Arj,A0) 


^/signal ' 

dAr\dA(ft ’ 


M(Arj, A<j>) 


Jljymixed 

dAr/dAip 


are the distributions for pairs from data and mixed events, respectively. The A (j) range is folded, 
i.e. for Atj) larger than n its value is recalculated as 2n — A(j>. In order to allow for a comparison 
with the RHIC and LHC results the pseudo-rapidity was calculated in the centre-of-mass (CMS) 
system. The transformation from the laboratory system to the CMS was performed assuming the 
pion mass for all produced particles. Electrons and positrons were removed by a cut on dE/dx the 
energy loss of the particle tracks in the TPC detectors. 


All charged, uncorrected 


All charged, correction 


All charged, corrected 





Figure 1: Example of the correction procedure for the correlation function of p+p interactions at 80 GeV/c 
beam momentum. Left plot shows correlation function C(At],A0) for uncorrected data, middle plot shows 
the binwise correction factors, right plot presents the corrected correlation function. For all plots the cut 
pr < 1.5 GeV/c was applied. 
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20 GeV/c, all 31 GeV/c, all 40 GeV/c, all 



Figure 2: Results on Ar)A0 correlations for inelastic p+p interactions. Results for all charged particle pairs. 
The correlation function is mirrored around (Arj.Atj)) = (0,0). 


3. Corrections 


In order to correct the measured raw correlation function for possible biases due to trigger and 
off-line event selection, track selection, etc. the same analysis was also performed on simulated 
data. The bin-by-bin correction Corr{Ar \, A(j )) was calculated as the ratio of the correlation func¬ 
tions for generated events from the EPOS [1] model (“pure”) and the same events after processing 
through Geant detector simulation and reconstruction (“rec”): 


Corr (A rj,A(j)) 


MCp Ure (A TJ , A(j) ) 
MC rec (Ari,A(j>) ’ 


(3.1) 


where MC pure is the correlation function obtained for generated events and MC rec is the correlation 
function for these events after detector simulation and reconstruction. For both “pure” and “rec” 
events the NA61/SHINE acceptance was applied. 

Corrected correlation functions (example in Fig. 1, right) were obtained by multiplying the 
uncorrected correlation function (Fig. 1, left) by the corresponding corrections (Fig. 1, middle), 
namely: 


C(Arj,A0) = C mw (Ar\,A<p) ■Corr(Arj 1 A(j)) 


(3.2) 


4. Results 

The corrected correlation functions for all charged pair combinations (unlike-sign pairs, posi¬ 
tively and negatively charged pairs) are presented in Figs. 2, 3, 4 and 5, respectively. 

Several structures can be seen in the plots: 
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20 GeV/c, + - 31 GeV/c, + - 40 GeV/c, + - 



Figure 3: Results on Ar\A(j) correlations for inelastic p+p interactions. Results for unlike-sign pairs. The 
correlation function is mirrored around (At],A0) = (0,0). 


20 GeV/c, + + 


31 GeV/c, + + 



80 GeV/c, + + 


158 GeV/c,+ + 


40 GeV/c, + + 






Figure 4: Results on A;]A@ correlations for inelastic p+p interactions. Results for positively charged pairs. 
The correlation function is mirrored around (Arj.Atp) = (0,0). 
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20 GeV/c, - - 31 GeV/c, - - 40 GeV/c, - - 



Figure 5: Results on Ar)A0 correlations for inelastic p+p interactions. Results for negatively charged pairs. 
The correlation function is mirrored around (Arj.Atj)) = (0,0). 

• A maximum at (Arj,A0) = (0.71), probably a result of resonance decays and momentum 
conservation. It is strongest for unlike-sign pairs and significantly weaker for same charge 
pairs. 

• Weak enhancement at (Arj,A0) = (0,0), likely due to Coulomb interactions (unlike-sign 
pairs) and quantum statistics (same charge pairs). 

4.1 Comparison with the EPOS model 

The ratio of the experimental results and predictions of the EPOS model is shown in Figs. 6, 7, 
8 and 9. In general EPOS reproduces results from p+p interactions well. However, the model does 
not reproduce the (Arj, Atj)) = (0,0) peak. This can be seen for positively and negatively charged 
pairs (Figs. 8 and 9 respectively). Note, that EPOS does not include Bose-Einstein correlations and 
Coulomb interactions. 


5 
















Two particle correlations from the energy scan with p+p interactions 


Bartosz Maksiak 


31 GeV/c, all 



A<|> [rad] 


A<]> [rad] 


A<[> [rad] 


80 GeV/c, all 


158 GeV/c, all 



A(j> [rad] 


A(|) [rad] 


Figure 6: The ratio of C(At],A 0) for data and EPOS. All charged particle pairs. The correlation function is 
mirrored around (A?],A0) = (0,0). 


20 GeV/c, + - 


31 GeV/c, + - 


40 GeV/c, + - 





Figure 7: The ratio of C(Ai ], A0) for data and EPOS. Unlike-sign pairs. The correlation function is mirrored 
around (Arj,A0) = (0,0). 
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20 GeV/c, + + 


31 GeV/c, + + 


40 GeV/c, + h 
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Figure 8 : The ratio of C(Ai7,A0) for data and EPOS. Positively charged pairs. The correlation function is 
mirrored around (Arj,A0) = (0,0). 


20 GeV/c, - - 


31 GeV/c, - - 


40 GeV/c, - - 



A(|) [rad] A<[> [rad] A<t> [rad] 


80 GeV/c,-- 158 GeV/c,-- 



A(j> [rad] A(|) [rad] 

Figure 9: The ratio of C(Arj.A(j)) for data and EPOS. Negatively charged pairs. The correlation function is 
mirrored around (Ar;,A0) = (0,0). 
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Figure 10: The effect of a transverse momentum cut pr <1.5 GeV/c on the measured correlation function. 


158 GeV/c, all charged, full p t spectrum 158 GeV/c, all charged, p y < 1.5 GeV/c 



Figure 11: The effect of a transverse momentum cut pj < 1.5 GeV/c on the correlation function obtained 
from the EPOS model. 

4.2 Influence of pr cut at 158GeV /c 

Originally, NA61/SHINE meant to study the effects of conservation laws, resonance decays, 
etc. on the correlation function ArjAtj). The effects of jets, though expected to be negligible, were 
reduced by the transverse momentum cut pr <1.5 GeV/c. An additional analysis was performed to 
check the effect of this cut. A comparison of the left and right plots in Figs. 10 and 11 demonstrates 
that the influence of the upper pj cut on the correlation functions for data and the EPOS model, 
respectively, are indeed negligible. 

4.3 Comparison with the ALICE experiment 

The NA61/SF1INE results are compared to the results from the ALICE experiment in Fig. 12. 
NA61/SITINE results show a stronger enhancement at A(j) ~ it and a significantly weaker maximum 
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at A(j) £3 0 (“jet peak”). 


20 GeV/c, all 


31 GeV/c, all 


40 GeV/c, all 




Figure 12: Correlation functions for all charged particle pairs from NA61/SHINE (five upper plots, same as 
in Fig. 2) compared with the ALICE results (three lower plots) taken from Ref. [2]. 

Figure 13 presents the energy dependence of C(0,0) (“jet peak”) and C(0. n) (“resonance 
hill”) from the NA61/SHINE and ALICE experiments. Monotonical behaviour is observed for 
both C(0,0) and C(0. n). C(0,0) rises with collision energy suggesting that the contribution from 
jets gets stronger with energy. On the other hand, C(0, n) decreases with collision energy. This can 
be explained by a decreasing contribution coming from resonance decays. 

Acknowledgements: This work was supported by the the National Science Centre, Poland grant 
2012/04/M/S T2/00816. 
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Figure 13: Correlation value C(Arj, A0) for “jet peak” and “resonance hill” as a function of energy in CMS 
frame. Left plot shows height of C(Ar) =0,A (f> =0) (“jet peak”). Right plot is a height of C(Ar) =0,A (f> =7i) 
(“resonance hill”). 
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